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Abstract 



The density of extended topological defects created during symmetry break- 
ing phase transitions depends on the ratio between the correlation length in 
the symmetric phase near T c and the winding length of the defects as deter- 
mined by the momentaneous effective action after a typical relaxation time. 
Conservation of winding number in numerical simulations requires a suitable 
embedding of the field variables and the appropriate geometrical implemen- 
tation of the winding density on the discrete lattice. We define a modified 
Kibble limit for the square lattice and obtain defect densities as functions 
of winding lengths in 0(2) and 0(3) models. The latter allows to observe 
formation of disoriented aligned domains within the easy plane. Their extent 
is severely limited by the momentaneous defect density during the course of 
the quench. 
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I. INTRODUCTION 



Formation of defects during symmetry breaking phase transitions has found increasing 
attention in a variety of applications ranging from condensed matter systems to cosmological 
scenarios [0 . The idea to consider baryons as topological defects in an iV-component chiral 
meson field <3>j constrained by $ 2 = const. 0, similarly has initiated attempts to estimate 
multiplicities of baryon and antibaryon production in high energy events or heavy ion colli- 
sions from the dynamics of defect formation during the cooling phase of an expanding hot 
hadronic gas Q. There are strong theoretical indications that during this cooling process 
the chiral O(N) symmetry is spontaneously broken near a critical temperature T c , and since 
the critical temperature is estimated to be of the order of the pion mass ( || and refs. 
therein), it may appear sufficient to describe this chiral transition in the framework of low- 
energy effective theory for the mesonic fields ||. However, perturbative methods for the 
evaluation of the relevant effective potential near the critical temperature are probably not 
very reliable, the role eventually played by other degrees of freedom is not really understood, 
so at present it is not even clear whether this chiral transition would be first or second order 
(see e.g. the discussion in ||). 

In analogy to phenomena observed in condensed matter systems there also have been 
speculations about the simultaneous ocurrence of extended domains with different average 
orientation of the aligned field ('disoriented chiral condensate' (DCC)) which rearrange into 
a uniform vacuum on a much longer time scale and therefore could manifest themselves 
through anomalous branching ratios for the production of differently charged mesons 0. 
Evidently, the existence of differently oriented domains and defect formation are intimately 
related, especially if the defects are of topological type which cannot simply disappear from 
a field configuration by local unwinding. 

Depending on the physical system under consideration the nature of the defects is de- 
termined by the choice of the manifold on which the fields live [10| : If the constraint $ 2 = 
const, is enforced, the fields live on the (N — l)-dimensional sphere S N_1 , conveniently 
parametrized by angular variables. This embedding can provide defects with a topologically 
protected winding number and thus in case of the chiral field allows to identify them with 
baryons. 

Near the phase transition, however, it may be appropriate to relax the constraint and 
allow all iV components of $ to move independently. Then a Euclidean R N embedding 
appears most convenient and in fact, the possible formation of DCC's in the chiral phase 
transition has been investigated in this framework ||TT|. The winding number of a field 
configuration embedded in R N is, however, not topologically conserved: it will undergo 
discrete changes, if at some space-time point all field components vanish, i.e. if the field 
configuration moves across the origin $ 2 (x) = 0. 

In order to preserve the identification of winding number with baryon number it is 
necessary to insist on the angular nature of the chiral field, while the additional degree of 
freedom is picked up by the modulus |$|, the scalar cr-field. In other words, the topologically 
trivial R embedding is replaced by the R l x S N ~ l manifold with nontrivial homotopy group 
Z. Here the origin $ = is excluded as a highly singular branch point where different angular 
sheets are tied together. Any field configuration which moves across this point, leaves a 
defect or antidefect of winding number ±1 at the spatial position where this happens, which 
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remains connected to the rest of the moving field by a string. The spatial structure of string 
and defect are determined by the dynamics of the classical field. The nontrivial structure of 
the modulus of $ along the string constitutes the 'bag', which interpolates from the vacuum 
value of the cx-field surrounding the defect to a (small) value in its center where the angular 
fields change rapidly from one sheet to the next. For a simple 1+1 dimensional 0(2) model 
this has been discussed in detail in |12| . 

These considerations naturally pose two questions: How do different embeddings affect 
the average multiplicities of defects and antidefects and the possible existence of disoriented 
domains during the chiral transition? And: What is the influence of the spatial structure 
of the topological defects (as determined by the temperature-dependent effective action) on 
their production multiplicities and on the possible formation of large extended disoriented 
domains? 

One may hope that for the essential features of defect formation the order of the transition 
is not very important. If the transition is first order then it would proceed through nucleation 
of bubbles ||13|| , which in their interior are characterized by an aligned chiral field with 



nonvanishing average value (|$|) = fo(T c ) of the chiral field as determined through the 
nontrivial minimum of the effective potential near T c . The orientation of the aligned field 
in different bubbles could be considered as random, so the average distance of their centers 
can be taken as an initial correlation length £(T C ) at T c . This length provides the relevant 
scale for defect formation, irrespective of its physical size. On this scale the total volume 
V considered can be identified with the total number of bubble seeds near T c . Growth and 
coalescence of bubbles then leads to the formation of topological defects within a typical 
formation time r characteristic for the damping of local fluctuations. 

Decisive for the multiplicity of the created defects is the ratio between their spatial extent, 
i.e. the 'radius' or 'winding length' l\y of the defects at the time of their formation and the 
correlation length £(T C ). The winding length is determined by the stabilization mechanism 
supplied through the effective action for the average field. 

If the winding length l\y is much smaller than the correlation length the defects can 



be considered as pointlike on the scale £(T C ). Then the purely combinatorial rules |14| of 
random lattices apply for their formation, i.e. their average density n is given by the Kibble 
limit 

n = {N + + N_)/V= l -(2)- D \ D (1) 

where (N + + N_) is the number of defects plus antidefects, D is the space dimension, and 
the factor A^ is the average ratio of the numbers of Z)-simplices and vertices in a large 
random lattice, i.e. Ai = l,Aa = 2, A3 = §§7r 2 , etc. |L5|]. In this case the random field 
fluctuations during the typical time r get transformed into a rather dense ensemble of defects 
and antidefects. These may become diluted in the following due to annihilation processes 
on a much larger timescale. 

If, however, the value of the winding length lw is larger than the average distance of the 
bubble seeds, then stabilization dynamics will prevent the formation of pointlike defects and 
antidefects inside the volume occupied by large extended stable defects or antidefects and 
(H|) has to be replaced by 
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n=(N + + N_)/V= 1 -(2 l -^)- D \ D . (2) 

i.e. before they have a chance to form, point defects and antidefects get instantly (i.e. within 
time r) eaten by big stable solitons which carry the net winding number of the chiral field 
within the volume occupied by them. The formation of these large extended solitons is a 
direct process which does not proceed through multiple annihilations of smaller defects. 

If for a sudden quench (£q <C r) the effective potential at T = is established before 
any dynamical or dissipative mechanism can affect the initial random field configuration 
then the winding length l\y{T) near T = will be relevant for the size of the 'cold' defects 
formed. With the relevant scale still given by the initial correlation length or average bubble 
distance £(T C ) near T c , the ratio lw(T = 0)/£(T c ) will enter into (@). 

If typical quench times are much larger than the formation time (tg ^> r) then the 
number of initially created defects and antidefects will be given by (Q) with winding length 
lw characteristic for T m T c . Subsequent changes in the defect multiplicities during the 
continuation of the quench depend on the scaling behaviour of their winding length: if 
lw{T) decreases with the decreasing temperature (as in three spatial dimensions where lw 
scales like /o(^) _1 ) then the (initially small) number of (initially large) defects stays small, 
so the multiplicity remains near the ratio lw{T c ) / £(T C ) in (0). On the other hand, if l\y(T) 
increases during the quench (as in one spatial dimension where lw scales like f (TY^ 2 ), the 
growing defects swallow many of the defect-antidefect pairs, so at the end of the slow quench 
the resulting multiplicities are the same as for a sudden quench, i.e. determined by the ratio 
l w (T = 0)/£(T e ). 

If the transition is second order no bubbles appear and the average distance between 
bubbles should be replaced by an instantaneous correlation length £(T). For finite rate of 
cooling £(T) remains finite near T c because critical slowing down prevents response over 



large distances [16]. On the other hand fo(T) stays close to zero which (in 3 dimensions) 
implies large lw- So also in this case the relevant winding lengths may exceed the pertinent 
correlation length of the plasma, so defect multiplicities can be much smaller than the 
limiting value given by eq. ([[]). 

Defect formation is intimately related with the size of DCC domains. If the transition 
proceeds in a way as to form a dense gas of pointlike defects and antidefects according 
to the Kibble limit (|IJ), then domains of topologically trivial mesonic field which fill the 
remaining space between the different defects consequently have to be small (in units of 
£). With increasing ratio lw/C the number of defects created decreases rapidly according 
to (0), but the chiral field in the spatial domains separating the defects is correlated due 
to the large winding length of the defects which 'organize ' the field surrounding them. In 
both limiting cases there seems to be little room for large (in units of £) domains of aligned 
but randomly oriented locally trivial field configurations. The actual physical size of these 
domains, of course, still depends on the magnitude of the basic length scale £ in physical 
units; this, naturally, depends on the physical system considered and on the order of the 
transition. For the hot chiral gas in the chirally symmetric phase above T c numbers for £ 
typically quoted range from 0.5 fm to 1 fm. On a lattice with a lattice constant of that size 
individual 'cold' baryons will cover just one lattice unit or less, while on formation near T c 
the 'melted' baryons [17j may smoothly extend over a large number of cells. We therefore 
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expect a situation where the above considerations about the relevance of the winding length 
typically apply. 

For a discussion of the chiral phase transition in the framework of effective meson fields 
it is therefore necessary to embed the fields in a manifold which allows for the topological 
protection of winding number, i.e. to use a separation into angular field variables and one 
modulus cr-field. Practically, in lattice simulations, where topological arguments no longer 
apply, conserving algorithms have to be used for the updating of configurations during their 
evolution. We will describe the results of such simulations for the ld-0(2) model in section 
II, for the 2d-0(3) model in section III. It is expected that the bias introduced into the 
initial configurations through explicit chiral symmetry breaking provides another efficient 
mechanism to suppress defect formation This effect should be discussed in connec- 

tion with the evolution of the a field. We shall concentrate here on the interplay between 
winding length and defect densities. Because the winding length is mainly determined by 
the dynamics of the angular field variables we postpone a detailed discussion of features 
related to the evolution of the modulus cr-field and, correspondingly, use unbiased initial 
configurations only. 



II. 0(2)-MODEL IN 1+1 DIMENSIONS 

Some of these features can be nicely visualized in the simple 0(2) model in 1+1 dimen- 
sions. As discussed above we concentrate here on the relation between winding length and 
defect density. To define the notation and for completeness we repeat some facts discussed 
in [O] concerning the topologically trivial R 2 vs. nontrivial R 1 x S 1 embedding. Defects 
arise as stable static solutions for the effective lagrangian density taken in the standard $ 4 
form (summation over a=l,2 and /i=0,l understood) 

£ ($) = \d^ a d^ a + ^($ 2 - ff - H&! . (3) 

The term H&i in (|3]) breaks chiral 0(2) symmetry explicitly. In order to keep the minimum 
of the potential for finite H at the if-independent value $ 2 = /q we define 

f 2 = fo-^r- (4) 

In accordance with conclusions from 1-loop renormalization group applied to (|3]) we will 
take A c and H as temperature independent, while fo(T) decreases as T rises from T = 
towards T = T c . 

The cartesian components $12 (—00 < $ 1)2 < 00) define the Euclidean R 2 embedding 
of this model, while modulus field a and the angular variable (0 < o < 00, — 00 < c/> < 00) 
related to <3>j through 

t) = a(x, t) cos0(x,t), ^ 2 (x,t) = a(x,t) sin0(x,t) (5) 

define the embedding into the S l x R 1 manifold. The a- and 7r-masses corresponding to 
these fields are 
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< = T> rnl = 2\ c fl + ml. (6) 
Jo 

A convenient unit for dimensionful quantities is the temperature-independent symmetry 
breaker H of mass dimension 2. So, in the following, if numbers are given for dimensionful 
quantities, they are to be multiplied by appropriate powers of H. 

In the Euclidean R 2 embedding the static part of ([5]) leads to stable nontopological 
solitons as long as the inequality 

/o^ > 21-28 (7) 

is satisfied. Otherwise, the only stable solution is the trivial ground state $i = /q, $2 = 0. 
This corresponds to approximately m a > 6.6 777/ 7j- . 3jS stability condition. For a typical 
coupling strength \ C /H = 60 (which we use in the following) this implies m n < 1.2 . On the 
other hand, in R 1 x S 1 , the winding number B of topologically nontrivial field configurations 
with fixed angular boundary conditions 0(±oo,t) = 27rn±, (n± : integer) 

1 f°° d 

B — — I —6(x, t)dx = n + — 7i-, (8) 
2ir J-00 ox 

is topologically conserved. For values of fo(T) below the limiting value (0) the angular field 
(f) of nontrivial static configurations collapses into pointlike defects at some position x = Xq, 
while the modulus field a describes a spatially extended 'bag' profile satisfying 

a" - \ c a{a 2 - f) + H = (9) 

with a(x — > x ) — > +0, and cr(\x\ — > 00) — > f . 

The radius of this bag scales approximately like ir/m^. Therefore, with / small near T c , 
where the stability condition (0) is not satisfied and the angular field collapsed to a point 
defect, the bag radius still is comparable to the angular winding length ly/ ps ir/m^ of the 
soliton, which solves the static 0(2) model (|3|), if the constraint $ 2 = fl is enforced, i.e., if 
we simply consider a U(l) model on a circle S 1 . 

On the other hand, if the stability condition (0) is satisfied, the static solutions of the 
0(2) model in R 2 and in R 1 x S 1 embedding coincide. Because the modulus field of these 
configurations stays close to f , the corresponding angular field essentially also coincides 
with the soliton of the model constrained to the circle S 1 with a = /o- 

Statistical ensembles of these stable solutions are expected at the end of a cooling process 
that evolves from an initial ensemble of random configurations with correlation length £. The 
quench is defined through the variation of fo{T) with time. For a numerical simulation the 
initial configurations before the onset of the quench we choose Gaussian deviates in the field 
components $ 12 (i), at the points Xi = i £, {i = 0, ..,N L ) of a spatial grid of total length 
L = Nl^ with lattice constant £ which defines the correlation length or distance of bubble 
seeds near T c . Periodic boundary conditions are imposed by choosing $1,2(0) = $i j2 (-^i,). 
From these random values of the cartesian field components the initial configuration of the 
angular field 4>(i) is obtained by taking the shortest path from the angle at point x^ to the 
angle at a neighbouring point x i+ %. This guarantees that the absolute value of the increment 
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in the angle from one point to the next is always less than ir. Due to the periodic boundary 
conditions in the cartesian field components the difference in the angular field 4>{Nl) — 0(0) 
then is an integer multiple of 2tt, which defines the winding number B of that particular 
random initial configuration. 

As we discussed above the winding length of nontrivial solutions in different embed- 
dings is similar to the winding length of the model with the modulus of $ constrained to 
<3> 2 = /o (T). The essentials of the relation between defect winding length and defect density 
therefore can also be observed within the simple U(l) (or nonlinear 0(2)) model. As an 
example we consider the sudden quench where the initially prepared configurations are ex- 
posed at time t — to the effective action with constant fo(T = 0). Then, within relaxation 
time r, the configurations evolve into an ensemble of stable defects with winding length 
l w « n/m n = 7t(/ /#) 1/2 . In fig.l the average final density n of defects plus antidefects 
averaged over an ensemble of 100 random initial configurations selected with total winding 
B = on an AT L =1000 lattice is plotted against the symmetry-breaking mass m n . 

0.25 
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0.15 
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0.05 
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FIG. 1. Defect plus antidefect density n = (N+ + N-)/Nl averaged over an ensemble of 100 
random initial configurations on an Nl = 1000 lattice after a sudden quench as function of inverse 
winding length for the nonlinear 0(2) model. The squares in the lower left corner show the 
same observable for the linear 0(2) model in i? 2 -embedding. 

The region where is less than tt is essentially characterized by a linear rise of the 
density n as expected from (|2]) for D — 1. Beyond this value the winding length becomes 
comparable to the lattice constant £ and we observe the onset of saturation with n slowly 
converging towards the Kibble limit (Jj) of n — 0.25 (for D — 1). A peculiar feature is that the 
increase in the number of finally surviving defects and antidefects with decreasing winding 
length proceeds in quite well-defined plateaus which are very robust against details of the 
evolution algorithm. Note that the absolute number of defects in the ensemble underlying 
fig.(l) is rather large: e.g. in the sudden decrease of n just below m n = 2.9 a total average 
number of 0.05 • 10 5 defects are eaten by a tiny increase of the winding length, while in 
the region (2.1 < m n < 2.8) mainly empty space in between defects is reduced. Evidently, 
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this is an indication of a lattice effect: as we have noted, if drops below it the defects 
start to grow beyond one lattice unit. As expected, for very small values of m n the average 
density is again rather smooth. Altogether, we may conclude that the experiment confirms 
the expectation expressed in (H). 

The squares connected by a full line in the left lower corner of fig.l show average defect 
densities as calculated in the unconstrained 0(2) model in Euclidean R 2 manifold. The 
theoretical limit of defect stability here (for A c = 60) is close to m n « 1.2. However, in the 
finite lattice calculation the defects disappear already near m n ~ 0.6. Beyond that value 
the resulting defect densities are zero, while in the region of stability the average densities 
closely follow the linear rise as obtained in the U(l) model. Irrespective of the stability 
condition, for the unconstrained 0(2) model embedded in R 1 x S 1 , the defect densities are 
similar to those of the U(l) model for all values of m n . However, in the region of instability 
where (^) is violated, average multiplicities in R 1 x S 1 embedding slightly exceed the U(l) 
results due to the smaller size of the radial bag (depending on the choice of A c ) (cf. eq. (|G|)). 

Relations between winding length lyy of defects and formation of DCC domains cannot be 
addressed within the 0(2) model, because the strength of the symmetry breaker necessary 
for finite lw prevents local alignment in random directions. Outside the range of the defects 
the angle <fi stabilizes always at multiples of 2ir. We therefore turn to the 0(3) model whose 
additional freedom allows to at least discuss partial aspects of DCC formation. 



III. NONLINEAR 0(3)-MODEL IN 2+1 DIMENSIONS 

Corresponding to our discussion of the constrained 0(2) model we consider in this section 
only that version of the 0(3) model where the modulus of $ is constrained to a fixed value 
/o. In two spatial dimensions /o is of mass dimension 1/2, so we use appropriate powers of 
/o as units for all dimensionful quantities. Then the fields $ are unit vectors $ which live 
on the sphere S 2 and conveniently are parametrized by two angles 9 and <fi 

$i (x, t) = cos </>(x, t) sin 6*(x, t) 
$ 2 ( x ; t) = sin 0(x, t) sin #(x, t) 

$ 3 (x,t) = cos#(x,t) (10) 

If $ converges towards a fixed unit vector, say e%, for spatial infinity on R 2 (sufficiently fast 
for the energy to converge) then the configuration space for smooth fields is disconnected due 
to the nontriviality of the second homotopy group ^(.S 2 ) = Z. Individual configurations 
then are characterized by the winding number B and local winding density p 

B = [ pd 2 x, (11) 



• (d^ x d,*) (12) 



(summation over spatial indices i, j understood). 
We consider the effective energy functional 
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m = J • + A 2 p 2 + — ($ - e 3 ) 2 ) d 2 x. (13) 

Evidently, A serves to scale the spatial extension of stable solutions, as it can be eliminated 
by x — > Ax. Both terms in S, the p 2 ('Skyrme') term and the explicitly symmetry-breaking 
last term, are necessary to stabilize solitons (the 'baby skyrmions' f[5|) with fixed winding 
length lw oc A. For our present purpose to study the relation between lw and average defect 
densities created in a symmetry breaking transition it is therefore sufficient to impose the 
quench through an appropriate time dependence of A. 

Discretizing the spatial coordinates to a two-dimensional lattice, the limit A C 1 will 
produce defects which are pointlike on the scale of the lattice constant, while A ^> 1 will 
create smooth density distributions which extend over many lattice units. A characteristic 
problem arises in numerical simulations on a discrete lattice if the spatial extent of solitons 
is of the order of a few lattice constants or even less: this implies, that differences in the 
angles 9 and <fi between neighbouring lattice points may be of the order of it. In that case 
the definition ( |12|) for the winding density p no longer applies. It is only for infinitesimally 
small differentials that the surface element dfl on S 2 which is the image of the elementary 
lattice cell d 2 x is given by dfl = Aarp d 2 x. But it is essential to maintain also for finite 
lattice constants the geometrical meaning of dQ as the surface area which the image of the 
unit lattice cell cuts out from S 2 . This can be readily implemented by defining 2np as 
the (oriented) area of the spherical triangle which is cut out on S 2 by the three geodesies 
which connect the end points of three <I>-vectors attached to three corners of an elementary 
lattice cell. This works for arbitrary relative orientations of these $- vectors and allows even 
to detect one complete soliton inside one elementary lattice cellQ. Clearly, this definition 
reduces to ([12]) in the continuum limit. Periodic boundary conditions (which compactify R 2 
to a torus), or $ — > e 3 (which compactifies R 2 to S 2 ) then lead to integer values of the total 
winding number B obtained by summing up all oriented spherical triangles. 

The topological considerations which guarantee conservation of B for continuous fields 
do not hold on the discrete lattice, because $- vectors attached to neighbouring points on the 
lattice can differ arbitrarily. However, 5-conservation can easily be re-implemented into the 
evolution of the configurations by allowing in each updating step only configurations^ which 
conserve B. In fact, this is a convenient way to compare evolutions which conserve B with 
others that allow for local unwinding of defects. So it is also not necessary to introduce spe- 



cific types of potentials (see e.g. ||20fl ) into (|T3|) to avoid the 'exceptional' configurations |2T 
as doorways for unwinding. 

For the first term in ( |13"D there is no similar compulsory extension to finite lattice con- 



In other words, we use a triangulation of the plane in terms of the triangles formed by two 
adjoining lattice links and one diagonal in each square lattice cell. Depending on the choice of the 
diagonal in the unit cell the area covered by the image of that cell on S 2 may differ by ±47r, i.e. 
the winding number may differ by ±1. However, for the ensemble average, this is not relevant. 

2 Practically, in a local updating procedure this requires to control the sum of p only in the local 
vicinity of the points which are affected by the update. 
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stants because it is not of geometrical nature. So, as usual we just interprete it as nearest 
neighbour interaction. 

Having thus defined the implementation of the model ( Jl3| ) on a discrete lattice which 
conserves those features essential for our present purpose, it remains to specify the initial 
configurations before the onset of the quench. Again we postpone here the question of any 
possible bias which may exist in these initial sets due to explicit symmetry breaking. Instead 
we choose a uniformly random set of ^-vectors on all interior lattice points. Only on the 
boundaries we impose $ = £3, i.e. 6 = with random cf). 

In the following sense this initial random set realizes the maximal average number of 
defects which can be accommodated on this square lattice: First it should be noted that 
by definition of our map (which maps a triangle formed by two adjoining orthogonal lattice 
links and the diagonal connecting their endpoints onto the interior of the smaller of the two 
complementary spherical triangles formed by the corresponding geodesies on S 2 ) the image 
triangle covers always less than half of the sphere, and thus by definition does not contain 
a defect. So, this is different from the Kibble mechanism which allows for defects inside 
elementary triangles by considering an additional lattice point inside the triangle. However, 
if we add the (oriented) areas covered by the two adjoining images of both triangles which 
make up the elementary square lattice cell then it may happen that this sum covers more 
than half of the sphere, and by definition we then may count that as one defect inside this 
elementary lattice cell. In fact, if A is very small, such a configuration will rapidly evolve 
into a configuration where the sum of these two spherical triangles (which represent the 
image of the lattice cell) is close to 4tt. We may call this the Kibble limit for the square 
lattice. It differs from the triangular limit ([!]) by a factor 1/4 (in two dimensions), which 
we interprete as a difference in the definition of the correlation length £ by a factor of two 
for uniformly random configurations on the vertices of a square lattice, as compared to the 
original Kibble counting which uses one additional point with random field asignment inside 
each of the two triangles which constitute the elementary square lattice cell. We consider 
the identification of £ with the lattice constant for a random asignment of field vectors to 
each vertex of the square lattice as the natural definition of the correlation length. 

This basic setup opens a wide variety of dynamical situations for systematic investigation. 
As just one example we shall discuss here only the case of a linear quench in the scale 
parameter A 

A = Xi + t (0<t<t Q ) (14) 

with Aj ^> 1 and Aj < 1. Starting from the initial random set, at the beginning of the 
quench a large value of Aj attempts to establish long range order with very small values of 
the local winding density p. When A approaches values aound A ~ 3 to 1, individual defects 
start to be formed out of this more or less organized low-density soup. Near the end of such 
a quench, for very small A we expext an ensemble of N + defects and iV_ antidefects each 
essentially localized on a single lattice cell. For t > tq the evolution can still be continued 
with A = A/ to allow for final stabilization of the defects formed, or for possible mutual 
annihilations on a much longer time scale. 

After the onset of the quench ordering of the configurations proceeds within typical 
relaxation times r and the extent to which the ordering can spread before the local defects 
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appear is governed by the ratio r/tQ, which is the appropriate measure for the quench 
velocity. For a fast quench we expect the final defect density to be close to our modified 
Kibble limit discussed above, because most of the initial random (incomplete) winding will 
directly evolve into complete local An winding. For slow quenches the random winding will 
have time to be smoothed away into large areas of low winding density, before finally a few 
defects and antidefects reappear. 

This naturally provides the appropriate setting to investigate formation and size of DCC 
domains. There is, however, a severe limitation to this discussion, also within the model ({13D : 
The stabilization mechanism of the baby skyrmions relies on a strong symmetry breaker, 
the last term in (|13"1). This term prevents formation of large areas of aligned but randomly 
oriented field configurations if we consider the complete field vectors In the spatial 
domains outside the range of the topological defects the field will essentially be aligned in 
£3 direction. Still, a partial aspect of DCC formation can be studied: In the so called 'easy 
plane' orthogonal to £3, the field components $1 and $2 can freely align in random directions 
and the partial correlation function 

C{\2- 1|) =< cos 0(1) cos 0(2) + sin 0(1) sin 0(2) > (15) 

(where the arguments 1,2 stand for two arbitrary lattice points) provides the information 
about the size of aligned domains. It should, however, be noted that the aligning force in the 
easy plane due to the first term in (|T3"D (the nearest neighbour interaction) is active only in 
areas where sin# ^ 0. This means that the formation of DCC domains can only be studied 
during the early parts of the quench where the strength of the symmetry breaker is not yet 
sufficient to form small localized defects with strictly aligned (in £3 direction) domains in 
between. As soon as these local defects appear (around A < 3) the disoriented aligned areas 
in the easy plane begin to re-randomize due to thermal fluctuations. So, the formation of 
DCC domains is closely connected to small but non-vanishing winding density, which by 
itself is an interesting aspect. 




-0.2 1 1 1 1 1 1 1 

2 4 6 8 10 12 

FIG. 2. Correlation function © C(r)/C(0) for fast (t Q = 10) (dashed line) and slow 
(tQ = 200) (full line) quench near A rs 3. 
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FIG. 3. <3?-field projected into (and re-normalized in) the easy plane for (a) fast (tg 
(b) slow (iQ = 200) quench near A f» 3. 
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FIG. 4. Winding density p during (a) fast (tq = 10) and (b) slow (ig = 200) quench near A « 3. 

In the following we describe some results of numerical simulations on a 51x51 lattice 
with an ensemble of random initial configurations which, however, have been selected as to 
satisfy 5 = 0. Of course, other values of B can be selected to study a possible influence of 
the total winding number. The algorithm conserves B as discussed above. For illustration 
we present two intermediate steps in the evolution of one specific arbitrary configuration 
during a fast quench (tq = 10 timesteps) and a slow quench (ig = 200 timesteps). 

The first set shows the situation shortly before \(t) crosses the value of A ~ 3, i.e. just 
before localized defects begin to form. The corresponding correlation functions fll5|) in fig. 2 



indicate that the ordering in the easy plane has hardly progressed beyond one lattice unit 
during the fast quench, while it extends beyond five lattice units after the slow quench. 
This becomes very evident if we look (in figs. 3) at the field $ projected into the easy plane 
and normalized to unity in the easy plane (i.e. the components of the arrows in the plot 
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FIG. 5. Final winding density p after (a) fast (tq = 10) and (b) slow (tq = 200) quench at 
A = 0.1. The resulting (anti)defect numbers are N+ = N- = 46 and N+ = N- = 15 in cases (a) 
and (b), respectively. 

simply are cos0 and sin0). This provides an instant view of DCC domains and defects 
and it shows that for the fast quench a high density of nascent defects leaves no space for 
large aligned domains while during the slow quench a well-developed DCC-pattern emerges. 
The corresponding winding density distributions p display a rough surface (fig.4a) versus 
essentially smooth areas of low density with few emerging local spikes (fig. 4b). Finally, 
figs. 5 show the resulting ensemble of defects shortly after the end of the quench where A has 
reached its final value of A/ = 0.1 and the local defects have stabilized. Counting the final 
number of defects and antidefects which evolve from this particular initial configuration we 
find N + = N_ = 46 for the fast (t Q = 10) quench and N + = N_ = 15 for the slow (t Q = 200) 
quench. 

Evolving M = 50 randomly chosen (but selected for B = 0) initial configurations through 
the quench ([P|) leads to an ensemble of defects and antidefects for each quench time tq 
as shown in fig.6. The modified Kibble limit as discussed above for the square lattice is 
iV + + iV_ — > (50 x 50)/16 ~ 156. The boundary conditions (9 = 0) reduce this value to 
about 144 for this lattice size. For the sudden quench (tq = 0) the measured numbers 
approach a mean value of only about 133 due to a few instant annihilations. For increasing 
quench time the defect numbers drop to a mean value of about 32 near tq = 200 as indicated 
by the full line in fig.6 which connects the mean values for the different quench times. The 
root mean square deviation from the average values decreases from about A « 10 for the 
sudden quench to A « 5 near tq = 200. Of course, these details depend on the specific 
choice of the quench mechanism used here. For applications to specific physical situations 
the appropriate temperature dependence of the effective action and temporal structure of 
the quench has to implemented. 
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FIG. 6. Final defect plus antidefect multiplicities for different quench times ig; their average 
values are connected by the full line. 

IV. CONCLUSION 

The purpose of this note is to study the influence of the spatial extent of topological de- 
fects on their formation during a dynamical quench in a symmetry breaking phase transition. 
Apart from the combinatorial factors of the Kibble mechanism it is the ratio of the momen- 
taneous magnitude of the winding length lw (as determined by the temperature dependent 
effective classical lagrangian) to the correlation length £ near T c which is essential for the 
resulting multiplicities of defects and antidefects. This leads to a crucial role of the quench 
time on the scale of the relaxation time r typically needed for defect formation. If during 
the early part of the quench the winding length is large on the scale of £, the resulting defect 
multiplicities are much reduced as compared to the combinatorial Kibble limit for pointlike 
defects. 

As illustration we presented numerical simulations in the nonlinear ld-0(2) and 2d-0(3) 
models which are characterized by a topologically conserved winding number B, in view 
of applications where this number corresponds to a physically relevant observable. For the 
lattice simulation it is essential to implement the geometrical meaning of the winding density 
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for the discrete map. Then boundary conditions on angular variables still serve to enforce 
integer values for B, although topological arguments for ^-conservation no longer apply. 
However, by explicit control during the evolution 5-violating steps can be excluded. If 
conservation of B is an important feature of the physical system under discussion, like baryon 
number conservation in chiral meson field models, extensions to the linear versions of O(N) 
models have to include the additional degree of freedom in the form of a modulus variable 
in an S 1 ^^ 1 x R l embedding. For comparison we have briefly discussed the topologically 
trivial i? 2 -embedding of the linear 0(2) model. Unfortunately, in these low-dimensional 
models, the winding length of defects is closely tied to explicit symmetry breaking. This 
leads to a dramatic difference between Euclidean and angular embedding if the symmetry 
breaking is sufficiently strong to destabilize nontrivial configurations. The same feature also 
hampers the investigation of DCC domains in these low-dimensional models. It is only in 
the easy plane of the 0(3) model where partial aspects of randomly oriented alignment 
can be observed. It has been argued that the bias introduced into the ensemble of initial 
configurations (for T > T c ) by explicit symmetry breaking causes another drastic reduction 
of defect multiplicities. This interesting aspect has been omitted here, again for the same 
reason that the symmetry breakers simultaneously have to serve as stabilizers for the defects. 

Of course, the real challenge is the formation of 3d-0(4) skyrmions as a nonperturbative 
model of how baryons emerge in the cooling of the hot meson soup. In that model stabi- 
lization of defects and explicit symmetry breaking through small mesonic masses are almost 
unrelated. Therefore all the interesting features concerning winding length, formation of 
DCC domains, effects of the modulus field, baryon number conservation, bias in the initial 
ensemble, all can be studied quite independently with the same methods used here for baby 
skyrmions. 
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